Induced long-time correlations in a two-component lattice gas 
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The distinguishability of at least two species of particles in the classical lattice gas with no 
interactions except hard-core exclusion entails additional interparticle correlations. A nonlinear 
mixing flow appears and manifests itself most pronounced in the case of significant difference between 
mobilities of species. It may result in the induced correlations for the slow component mediated 
by the fast one. In the quasi-one-dimensional case, the long-time correlations are demonstrated 
to take place in the slow component, that is similar to the hydrodynamic correlations between 
colloidal particles. In the adiabatic approximation, these correlations may come into play only in 
O-i the non-equilibrium case with the flow of the fast component present in the system. 



I. INTRODUCTION 

Multicomponent systems often exhibit unusual behavior and transport effects Many of 

these effects are caused by the additional correlations which are associated with particle distinguish- 
ability. The short-range "back correlations" (memory effect) [l| in the classical two-component 
lattice gas with excluded volume constraint may serve as an example. At the macroscopic level the 
particle distinguishability manifests itself as the additional nonlinear flux which is responsible for 
the mixing of different species of a lattice gas. This flux leads to the array of anomalous diffusive 
transport effects, e.g., the drag effect, formation of the drifting spatial structures 0-131' effects of 
"negative" mass transport [1, Q . 

When the difference between the values of the component mobility is significant, the nonlinear 
flux can lead to indirect correlations in the slow subsystem. These correlations are induced by the 
light component that plays the role of interaction mediator for the particles of the slow subsystem. 
Such a phenomenon is similar to the hydrodynamic interactions between colloidal particles 
One of the intriguing effects observed in these suspensions is the presence of long-time correla- 
tions between liquid-embedded particles when the velocity correlation function of a particle pair 
is characterized by the long-time negative tail [Tol-[T2| . 

In this paper we are aimed to show the possibility of similar long-time correlations for the slow 
component of a lattice gas. For this purpose, we consider the simplest case of a two-component 
lattice gas by using the adiabatic and the uniform clusters approximations. 



II. CORRELATIONS IN A TWO-COMPONENT LATTICE GAS 

We consider the simplest two-component lattice gas model subject to an excluded volume con- 
straint. The kinetics is defined by the jumps of particles to the neighboring vacant sites. The varia- 
tion of the i-th site occupancy by the particles of sort a during the time interval At, tq ^ A< ^ r; 
(tq is the duration of a particle jump to a neighboring site, t; being the life time of a particle on a 
site), is described by the standard continuity equation (see, e.g., [ll[T3|) 

<(t + At) - <(t) = 5] ( - J^) + (1) 

j 

where nf = 0,1 are the local occupation numbers of particles a at the i-th site, J," = uf^n^hj/S.t 
gives the mean number of jumps (of particles a from site i to a neighboring site j per time At), 



^ij — Va is the frequency of these jumps, hj — 1 — X]/?'^^ the occupation number of the jth 
site vacancies. The term SJf stands for the Langevin source of the ith site occupation number 
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fluctuations. Equations for the average local occupation numbers may be obtained from Eqs. ([T]), 
using the local equilibrium approximation (Zubarev approach) [l^, which in our case coincides 
with the mean field approximation [l5| . In the case of two components, introducing time derivatives 
in Eqs. ([T]), see [l^, the macroscopic equations take the form 0, 0, 

i^m^m = V^m + V{mVn — nVm) + 5Jrm (2a) 



v-^n = V'^n-V{mVn-nVm) + 5Jn, (2b) 

where m and n are the average occupation numbers of the two components at the point r, Vm 
and Vn are the jump frequencies of particles m and n, correspondingly. These equations are 
obtained in the long-wavelength approximation assuming that the characteristic scale of the gas 
inhomogeneity is much grater than the lattice constant a. The dimensionless spatial coordinate 
r — > r/a is introduced. 

Distinguishability of the particles in the noninteracting lattice gas leads to the additional nonlin- 
ear flow nVm — mVn in the macroscopic Eqs. ^ that is absent in the case of a one-component gas. 
An important condition for the onset of this additional term is the excluded volume constraint, 
i.e., when a site may be occupied by one particle only. This flow is responsible for the mixing of 
different particle species and leads to additional correlations in the system. 

In the case of considerable difference between the species mobilities, Vn/i^m ^ 1, the indirect 
correlations can occur in the slow subsystem n. These correlations are induced by the fast com- 
ponent that plays the role of the "interaction mediator" . Such correlations are similar to the 
hydrodynamic correlations between colloidal particles [ll| and may have the long-time character. 

In order to show the presence of the long-time correlations in the heavy subsystem of the lattice 
gas we consider the fluctuations of the heavy particles drift caused by an external field. They 
manifest itself as the fluctuations of hydrodynamic velocity field u(r, t) . In what follows, we assume 
that the characteristic time-scale Tu of the fluctuations u(r,t) satisfies the condition Vm ^ ^ 
Vn and restrict ourselves to the adiabatic approximation, Vm -t- oo. It means that the distribution 
of the fast component instantly gets into the stationary one with the momentary distribution of 
the slow component. Then, kinetic equations take form 

= V^m -f V(mVn - nVm), (3a) 
Vn^n = V'^n -\/{mS/n~nVm) + Vn^u{r,T)Vn, (3b) 

The light component m(r, t) depends on a "slow" time t as on a parameter. It should be noted, 
that the fluctuations of drift are associated with the fluctuations of the particle jump probability 
I'Ti+i — + between the neighboring lattice sites i, i + 1 and in general case this term may 
have a more complicated form. For example, in the case of the asymmetric jumps of the lattice 
gas particles, v^i_^.l = Vn + u, ffj.^ — Vn, the drift term reads v~^u{\/n — Vn^ — \/nm). However, 
here we restrict ourselves to the case given by Eqs. ^ and, similarly to consider a quasi- 
one-dimensional situation when interparticle correlations become most significant. In this case, 
one can easily eliminate the fast variable m from Eqs. ([3]), see @, and obtain the self-consistent 
equation for the heavy component n 



dn{x, t) 
dr 



/ T/ \ ds \ d'^n(x,T) , .dn(x,T) 



J(r) = -(^,-^0 (^J^ ^,_X,rw) • (5) 

Here we introduced the dimensionless time t = i/nt. The quantity J(t) [Eq. ([5])] corresponds 
to the average macroscopic flow of the light component m{x,T) passing through the system 
m{—L, r) = TOi = const, m(L, r) = 7712 — const, and n{—L, t) = n{L, r) = being the values of the 
mean occupation numbers m and n at the system boundaries x — ^L. 

It follows from Eq. @ that the correlations in the slow subsystem appear only in the nonequilib- 
rium case with the nonzero flow J(t) that is responsible for the nonlinear and nonlocal character 
of the diffusion process. 
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Next we make use of the uniform density cluster approximation that is often exploited, e.g., for 
the description of Coulomb explosions of deuterium clusters [itI , [isj . Under this approximation, 
the density distribution n{x, r) is considered as a system of uniform one-dimensional clusters. It 
is justified when the cluster size is much smaller than that of the system. Each such cluster is 
characterized by its "center of mass" rfc(r) and width Rkir), so that the distribution of the slow 
component can be written as 

^^-f Mr) = Qk/Rkir) if \x - rfc(T)| < i?fc(T)/2, 

' ^"\0 ii\x-rk{T)\>Rk{T)/2, W 

where Qk is the number of the particles contained by the fcth cluster. Further, we assume the 
characteristic spatial scale ^ of the fluctuations u{x, r) to be much grater than the typical cluster 
size but less than the distance between them, Rk ^ ^ <C \rk+i —rk\. Then, evaluating the moments 
of distribution ([B]) from Eq. ^ one can obtain the equations for rfc(r) and Rk{T) 

= J(t)GK(r)] - :^-Wir), (7a) 
^ = 2D[rk{r)] + J{r) ( ^ QiG[m{r)] - ^ QiG[m{T)]\ . (7b) 

\l<k l>k / 

The first equation ((7a|) describes the drift velocity fk of the fcth cluster. The motion of the 
cluster as a whole is induced by the fluctuations of the external field Uk and by the fiow J(r) = 
2Lj{2L + J2kQkG[nk{T)]}~^ of the fast component (the drag effect). The function G[rifc(r)] — 
[1 — nfe(T)]~^ + [1 ~ 'n'k{T)]~'^ describes the drag velocity slowdown caused by the cluster expansion. 
The second equation ((7b| describes the rate of the cluster expansion that is defined by the diffusion 
coefficient D[rk{T)] = 1 — (mi + m2)/2 + J{T)rk{T) to be cluster position dependent and by the 
drag velocities JG{nk) of other clusters, see ([7a|) . Note that the quantity j — — Am/ {2L) ^ 1 is a 
natural small parameter since |Am| — |m2— < 1. 

In fact Eqs. Q describe the Brownian-like dynamics of the penetrable, expanding clusters em- 
bedded into the flow of the fast component. For simplicity, next we consider the system of the 
two initially identical clusters, Q — Qi{2): R{^) — ^1(2) (0)i with the ^-correlated fluctuations, 
{uuIt)ui{t')) =T5m5{t-t'). 

We are interested in the response of one cluster caused by the initial velocity perturbation of 
another cluster. For this purpose, we consider the pair correlation function of velocity fluctuations 
Ci2(t) = {5vi{t)5v2{Q)) = ([ri(r) - wi(r)] [7-2(0) -^2(0)]), Ufe(T) = rfc(r), where ffc(r) satisfles 
unperturbed Eqs. ([7]) (i.e., Eqs. ([7]) provided u\^{t) =0). Supposing that the fluctuations Mfe(t) 
are small, v~^\u\^\ <IC |ufc(i)|, and using the perturbation theory in v~^\u]^\ and j (we assume, that 

< v~'^\uk \ < i), one can get 

-2.3 24rQ / R{t) 2QR{r) GoRjr) \ 

Oi2irj^^„j I?o U^W-Qr [^M-QP [i?(r)-QP;' 

where i?(r) = [24L>t -I- i?2(0)]i/2^ Do ^ I ~ (toi +TO2)/2 and Co = [i?(0) -h g]/[i?(0) - Q]^. 

The correlation function ([5]), Fig. (P), has the long-time negative tail with the asymptotic be- 
havior Ci2(t — j> 00) ~ —T~^. The long-time correlations between the two diffusive clusters are 
similar to hydrodynamic ones appearing between colloidal particles confined in a liquid-filled linear 
channel 



III. SUMMARY AND DISCUSSION 



The significant difference of particle mobilities in a two-component lattice gas leads to the induced 
correlations in the slow subsystem mediated by the fast one. They manifest itself as a negative 
long-time tail in the behavior of the pair correlation function and are similar to the well-known long- 
time correlations between colloidal particles [ill, [13 • Note that there is no additional interparticle 
interaction introduced in the lattice gas except that due to hard-core exclusion. In order to 
demonstrate this phenomenon we have considered the one-dimensional case and used the simplest. 
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FIG. 1: Velocity cross-correlation function of a cluster pair as a function of time [Q — 0.5, -R(O) — 1, 
Am = 1, 7 = 24rQj^/(Do!^^)]. The inset shows the long-time asymptotic behavior, the dashed line 
indicates a slope of 1. 



somewhat rough, approximations. In particular, we have used the adiabatic one and neglected 
the retardation effects. As a result, the correlations instantly spread all over the system and the 
averaged flow of the fast component does not depend on the distance between the clusters of 
the slow one. It give rise to the weak dependence of the correlation function on the inter-cluster 
distance that make itself evident in the higher orders of the perturbation theory expansion. 

Note that the correlations in two- and three-dimensional cases may be considered in the same 
manner by using the uniform density cluster approximation. 

Considered correlations may be associated with the "elastic-type interactions", similar to the 
case of the liquid-embedded particles [ll,!!^]. However, another type of the induced correlations in 
the slow component may be caused by the equilibrium fluctuations in the fast one and associated 
with so-called Casimir-like forces f2l)], that is beyond the scope of this paper. 
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